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Abstract 

Integrable problem of three vortices on a plane and sphere are considered. 
The classification of Poisson structures is carried out. We accomplish the 
bifurcational analysis using the variables introduced in previous part of the 
work. 



1 Introduction 

In this part of our work the integrable problem of dynamic of vortices on planes 
and sphere are considered. We base on Hamiltonian form of reduced systems, which 
poisson structure was given in the previous part [TH] . 

The problem of motion of two vortices on plane was completely investigated by 
G.Helmholz [T3|, who had established, that in general case two vortices make an 
uniform rotary motion around of vorticity centre, with frequency 

2vrr 2 ' 

where r*. is intensity of point vortices. The vorticity centre thus is moving uniformly 
and rectilinearly. If T\ = —T2 the vorticity centre is situated on infinity, and two 
vortices are moving forward. 
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The problem about a motion of three vortices is much more complex. Despite 
of numerous works, first of which were the dissertation of Grobli 1877 [2] and 
Greenhill's research jS], the complete and evident classification of a motion still 
does not exist. From the point of view of integrability problem, this problem was 
considered by Poincare in the treatise on the vortices theory He wrote out 
explicitly a complete set of noncommutative integrals. In a modern period the 
problem of a motion of three vortices on a plane was studied in works [121 E3 E3 EE 
|2*U] , and from the point of view of the topological analysis — in . Unfortunately, 
these works have added a little to achievements of classicists both in presentation, 
and in completeness of the description of motions. Their basic contents come down 
either to geometrical interpretation of a motion and computer modeling of separate 
trajectories, or to some general topological constractions, that is not coordinated to 
physical behaviour of system. Partly, it is caused by the fact, that the problem of 
three vortices on a plane does not belong to that integrable systems, which complete 
analysis is possible in a class of rather simple (for example, elliptic) special functions 
(the general solution has indefinitely-sheeted branching on complex plane of time 
because of logarithmic terms in Hamiltonian) . The exception consists of some special 
cases (for example, the case of equal intensities of vortices). In comparison with a 
problem of three vortices on a plane, a motion of three vortices on sphere, that is 
also being integrable, is not investigated at all. 

We give here new analysis, based on algebraic and geometrical research of the 
reduced system (in variable, which in the previous part of work ^H] were named 
"internal", as though we are not adhere this term anymore), and then we analyze 
absolute motion, without use of exact quadrature. Such approach, basing on rep- 
resentation of the equations of motions on algebra 18], allows to establish some 
analogies between problem of three vortices and Euler-Poinsot case in rigid body 
dynamics, and also to receive more evident description of motions of system. 



2 Algebraic classification 

As it was shown in the previous part of work ^H] , the equation of a motion of three 
vortices on the plane can be written down as Hamiltonian system, determined by 
Lie-Poisson brackets in form of 

{Mi, Mj} = -4a fe A, 

{Mi, A} = (a, - a k )Mi + (a,- + a fc )(M, - M k ), 
and by Hamiltonian ^ 

h = -— (r 2 r 3 iogM 1 + r 1 r 3 iogM 2 + r 1 r 2 iogM 3 ). (2) 

87T 

where M k are squares of pairwise distances between vortices, A is a size of the 
oriented area pulled on three vortices (here and further we shall assume, that the 



2 



indexes i,j, k accept accordingly values 1,2,3 and their cyclic rearrangements, and 
the coefficients 01,02,03 are inverse intensities a k = 1/T k ). 

Lie-Poisson brackets (0) are degenerated and have two central functions. One of 
them (linear) is integral of the complete moment 

D = J2 a kM k , (3) 

k 

the other (square-law Casimir function) arises from geometrical Henon ratio, that 
connects the area of triangle with sides 

F = (2A) 2 + M 1 M 1 + M 2 M 2 + M 3 M 3 - 2(M 1 M 2 + M X M 3 + M 2 M 3 ). (4) 

For real motions F — 0. 

The real type of Lie-Poisson algebra (JTJ depends on values of intensities Ti, T 2 , r 3 . 
In fact, choosing a new basis in form D, ei, e 2 , e 3 , where D is defined by © 



ei 



(a 2 - a 3 )Mi + (o 3 - ai)M 2 + (01 - a 2 )M 3 

62 = vsi • (5) 

(a 2 a 3 - a 2 )Mi + (a x a 3 - a 2 )M 2 + (aia 2 - a 3 )M 3 
e3 — 
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where 

v4 = |aia 2 + a 2 a 3 + aia 3 |, 
B = (ai — a 2 ) 2 + (a 2 — a 3 ) 2 + (ai — a 3 ) 2 , 
under the condition + G2G3 + aiQ3 ) > 0> ( 6 ) 

we receive that the algebra of vortices is decomposed in the direct sum 1(4) ps 
R © so(3) : 

{£>, e fe } = 0, {e x , e 2 } = e 3 , {e 2 , e 3 } = ei {e 3 , ej = e 2 , (7) 
and by condition ( a%a2 + a2fl3 + aiQ3 ) < , (8) 

in direct sum 1(4) ^ IR © so(2, 1) : 

{D,e k } = 0, {ei,e 2 } = -e 3 , {e 2 ,e 3 } = ei, {e 3 ,ei} = e 2 . (9) 

Though for case equal intensities the coefficient B = 0, using the passage to the 
limit, it is easy to show that basis (0) can be correctly determined in this case as 
well. 



3 



Let's consider, what features of a motion can be obtained from the structure of 
such algebraic decomposition. Using squared Casimir function for algebra so(3) 



G = e{ + e 2 + ei (10) 

and expressing in it a vector e\ with the help of Heron ratio (@J) with F = 0, we 
discover, that relative dynamics of vortices with intensities, satisfying to a condition 
(JBJ) will be equivalent to a motion of some "representing" point on symplectic sheet, 
that is a two-dimensional sphere with radius, determined by value of linear Casimir 
function 

2 2 2 1 i Oi0 2 (2 3 _D 



e( + ei + ei = — [ . (11) 

16 V a x a 2 + a 2 a 3 + a x a 3 ) 

The situation for algebra so(2, 1) is similar that arise in a case (jHJ). Using Casimir 
function of algebra so(2, 1) 

G 2 = el-el-el (12) 

it is possible to notice, that the relative motion of vortices is reduced to the motion 
of point on surface of two-cavities hyperboloid, that is determined by condition, 



1 / ctia 2 a 3 .D 



16 V a\d2 + a 2 a 3 + (I1CL3 



(13) 



In case of (fU), motions of representing point (and therefore three vortices in system 
of vorticity centre) is finite for every D, therefore this case we shall futher name 
"compact", and in case (JHJ), with which obviously may be exist running up trajec- 
tories by additional conditions for Tj and for initial vortices positions, the name 
"uncompact" . 

Motions, arising under condition 

(ai<2 2 + a 2 a 3 + aia 3 ) = (14) 

require separate consideration. In basis 

_ aiMi + a 2 M 2 + a 3 M 3 
(01 + a 2 ) 2 

A 

ei 



(ai + a 2 

e 2 



(15) 



2(ai + a 2 ) ' 
(ai + a 2 + a 3 )M 3 - (Mi + M 2 )(ai + a 2 ) 
63 — 



2( ai + a 2 ) 2 

Lie-Poisson brackets can be reduced to form 

{D, e k } = 0, {e x , e 2 } = e 3 , {e 3 , ej = — D, {e 2 , e 3 } = ei. (16) 
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In this case algebra of vortices is not decomposed in the straight sum of one- 
dimensional and three-dimensional algebras, and becomes four-dimensional solve- 
able algebra with maximal solveable ideal N = {D,ei,e 3 }. Algebra has squared 
Casimir function, following from (j3J) 



Real dynamics of vortices occurs on symplectic sheet, that is set by conditions G = 
and D = const, which determine paraboloid, taking place through the beginning 
of coordinates. By D — 0, being necessary collaps condition (merge of vortices), 
paraboloid degenerates to the straight line, conterminous with an axis e 3 . 



3 Bifurcational analysis of motion of vortices on 



Let's describe most evident geometrical interpretation of motions, used in [TH] and 
presented on Fig. 1. In space Mi, M 2 , M 3 the level of linear integral of the moment 
sets a plane. The inequalities > allocate on it the area, in which motion 
occurs. It is necessary to exclude from this area unphysical values of distances, for 
which an inequality of a triangle (A 2 < 0) does not hold. This area is shown at a 
figure by black colour. When approaching to it, in the equations of a motion for 
variable Mi, M 2 , M 3 , we schould change the sign of time and direction of motion. It 
is necessary to note, that the described geometrical interpretation is not equivalent 
to a motion of representing point on symplectic sheet, and changes of direction of 
motions are the consequence of features of a projection from space A, Mi, M 2 , M 3 
in space M 1 ,M 2 ,M 3 . 

Let's hold condition then bracket (£[]) is defines algebra so(3)©R. In this case 
symplectic sheet is compact (§ 2 ), and motions of vortices are finite. Let's assume, at 
first, that signs of all intensities are identical — in this case condition © is obviously 
fulfilled. The first integrals of a motion of vortices on a plane we shall write down 
as 



D = aM + a 2 M 2 + a 3 M 3 , 

(18) 

h = Mt 1 M?M 3 a \ 

(integral of energy, for convenience, is presented in exponential form). The integrals 
(I18|) are dependent, i.e. Jacobi matrix of the first integrals (I18j) is degenerate 



only in one case: Mi = M 2 = M 3 (in three-dimensional space .\/ ; . M>. M :i the 
vortices will form correct triangle). Curves, resulting with crossing of levels of 



G 2 = e\ + e\ - 2De 2 . 



(17) 



plane 




(19) 
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integrals (|18|) are similar to polhodes in rigid body dynamics (Poinsot interpretation 
of Euler case) and, in a compact CcLSG, £1X6 ovals (Fig. 1). The solution as correct 
triangle occurs with a contact of a surface h and plane D in one point. They limit 
by energy the area of possible motions (APM) from above and the corresponding 
bifurcational curve has a form ,^ „ n 

H(D) = fc) ' (20) 

Remark 1. The particular solutions, appropriate to the given curve — three vor- 
tices in tops of correct triangle, rotating as a rigid body around of vorticity centre 
are refered as "Thomson's" and are steady with fulfilment of a condition (Q)J). J. 
J. Thomson has specified them for any number of vortices of equal intensity, and 
has shown, that in linear approach such configurations will be steady for number 
of vortices N ^ 6, and with N ^ 7 — unstable (Thomson theorem) However 
because of presence of resonances in system the linear approach is not sufficient. 
The analysis of stability in nonlinear approach with use of Birkhoff normalization 
was carried out in [6]. It was found, that the theorem Thomson is fair in the exact 
statement — in sense of Laypunov stability. 

The other change of motion occurs with a contact of a crossing line of integrals' 
levels with curves, determined by the equation A (Mi, M 2 , M 3 ) = (y/Mi + \/ Mj — 
v^Mfc), that arise from inequalities of triangles and limit physically allowable area 
of values (see Fig. 1). 

Points of a contact correspond to collinear configuration of three vortices. Three 
vortices thus settle down on one straight line and rotate as a unit around of vorticity 
centre. 

Remark 2. Collinear and the triangular configurations in dynamics of three vor- 
tices have analogues in the classical selestial mechanics \1 If . There are, accordingly, 
Euler and Lagrange particular solutions of three body problem. 

Collinear configurations are obtained from a condition of a contact 

A = 0, A = 0. (21) 

that lead, having same sedate functional dependence, which in a general case can 
be presented as 

H(D) = / (ai, a 2 , a 3 )D^ +a * +a >\ (22) 

where f(a±, a 2 , a 3 ) some function from parameters. Passing to homogeneous coordi- 
nates x, y 

Mj •'•.!/,, M 2 = yM 3 , y =(l + x)\ x = z\ (23) 

we discover, that collinear configurations correspond to positive roots of the equa- 
tions of the third degree 

ai(l ± z)(z ± 2) + a 2 ((l =F z)z - a 3 (l ± 2z)z{\ ± z) = 0, (24) 
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where 

ai = r 2 + r 3 , a 2 = ri + r 3 , a 3 = ri + r 2 . 

The lines of crossing of levels D and h can concern borders of area of physically 
possible motion A = in space of variables (Mi,M 2 ,M 3 ) by various ways (Fig. 1). 
Depending on it the quantity of bifurcational curves also will be different. In partic- 
ular, in a case, when all of intensities are different a\ 7^ a 2 7^ a 3 7^ a\ three branches 
corresponding to collinear motions will exist (Fig. 2a). With equality of two of them, 
for example, a\ = a 2 7^ a 3 two appropriate bifurcational curves merge. And, at last, 
in case of equal intensities CL\ — (2 2 — Ct.3 all three curves degenerate in one. 

The calculation of angular speed of Thomson's and collinear configurations con- 
cerning vorticity centre is given, for example in ^3] . It monotonously decreases with 
increase of the complete moment of system vortices. 

As a measure of stability of stationary configurations it is possible to use a square 
of untrivial eigenvalues A 2 of linearized system of the equations of a motion of three 
vortices. For Thomson's configurations it is simple to receive the explicit formula 

,2 o (a 1 + a 2 + a 3 ) 3 

A = -3aia 2 a 3 — , (25) 

that shows, that in a considered case they are steady, in difference from collinear, 
which as shown on Fig. 2a a are unstable. 

Generally speaking, the following statement is fair, that is confirmed by geomet- 
rical interpretation on Fig. 1: If Thomson's solution is steady, the appropriate set 
of collinear solutions are unstable, and on the contrary. Therefore condition (|25|) 
determines a type of stability not only Thomson's configurations, but also collinear 
configurations with the given values of D. 

Let's consider other case of a motion of three vortices, with fulfilled conditions 
© (easy to see, that all other possible cases, when the condition is fair, could be 
reduced to two considered cases). Let's assume, that the intensity of one of vortices 
has opposite sign in comparison with other two (for example, Ti = 1/ai < 0). The 
condition © in this case means, that — Fx < T 2 + T 3 (i.e. intensity of chosen vortex 
is greater than the intensity of other two). 

Let's specify only main differences of this problem from previous. First of all 
from the formula (J25j) for Thomson's configurations it follows, that they are unstable. 
Collinear configuration (it is not three of them anymore, but only one) thus, on the 
contrary, is steady. Bifurcational curve of appropriate unique collinear configuration 
(see Fig. 2b) is above a curve, appropriate to Thomson's configuration (it is simple 
to show, that under condition of © complete moment is always more than zero), 
and the area of possible motion coincides with the whole square D > 0, h > 0. 

To obtain more complete representation about the considered motions, let's 
study a phase portrait of system of three vortices on a plane in symplectic coor- 
dinates. 
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4 Symplectic coordinates for vortices on a plane 



Since symplectic sheets represent two-dimensional spheres f|18|) . it is convenient to 
use cylindrical coordinates 1,1 as canonical coordinates (we shall keep for them a 
designation of Andoyer-Deprit coordinates, adopted in rigid body dynamics jTHj) 
for algebra so (3) : 

ei = L 

e 2 = V G 2 -L 2 sinl, (26) 
e 3 = \JG 2 — L 2 cosZ. 

Expressions for variables Mi, M 2 , M 3 in canonical coordinates look like 

J - (27) 



where 



aj(aj - Oi) + a k (a k - Oi))\/ r cos/ 



A = ai<22 + 0,2(33 + 0-304, 

5 = (ai — a 2 ) 2 + (a 2 - a 3 ) 2 + (ai - a 3 ) 2 , 



F = VG 2 - L 2 . 

In that special case of equal intensities the derived expressions (|T2"j) lead to the 
expression for Hamiltonian mentioned in ^21 E] In last work it is used the classical 
way of a reduction to a problem with one degree of freedom for arbitrary intensities 
(such as Jacobi's unit elimination). 

The phase portrait in canonical coordinates describes a motion of representing 
point (L, I) on a surface of two-dimensional sphere. The pulse L has a value of 
the oriental area parallelogram, that is constructed on three vortices (L = 2A). 
Development of a phase portrait is presented on Fig. 4a-d for four different specific 
cases of the intensities of vortices ratio. 

The given figures illustrate the analogy between the motion of three vortices 
and rigid body dynamics. Comparing Fig. 2a (in case of equal intensities) with 
phase portrait of a Euler-Poinsot problem (see for example JE]), it is possible to 
connect a collinear configuration (situated on a straight line L = 0) with unstable 
permanent motions of a rigid body around of an middle axis of ellipsoid of inertia, 
the Thomson's solutions (with L/G = 1) — with rotations around of a large (small) 
axis of ellipsoid of inertia. The periodic solutions of a problem of two vortices (two 
of three vortices always are in one point, and their intensities are added), situated on 
straight line L = it is possible to connect with steady permanent rotations around 
of a small (large) axis of elipsoid of inertia. When system pass through the collinear 
state ( three vortices on one straight line — L = 0), oriented area changes the sign. 
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Deformation of a phase picture with unequal intensities is shown on Fig. 4b, c. Thus 
Thomson's solutions are displaced with straight L/G = 1, that on sphere should be 
imagined as being pulled together in a point). Fig. 4d present the phase portrait, 
with one negative intensity, but condition (jUJ) is also fulfilled. 

Other motions of system of three vortices, that is not a state of equilibrium on 
phase portraits, correspond to quasiperiodic motions of vortices on a plane. Depend- 
ing on the topology of a phase curve these motions have various features — some of 
them repeatedly pass through the collinear state, others — never get into it. Such 
areas are separated by separatrix curve, for that curve three vortices are approach 
to collinear (unstable) configuration infinitely long. The particular trajectories of 
vortices are given in various works [HH IT51 IT^ ITU] , beginning from Grobli j2] and 
Greenhill jS]. 

Appropriate simplectic coordinates for algebra so(2, 1) could be chosen as 

ei = L 

e 2 = VG 2 + H shl, (28) 
e 3 = VG 2 + L 2 chl. 

where G, as well as in a case so(3), parametrize various simplectic sheets. The dif- 
ference of the case (JEJ is displayed in noncompactness of symplectic sheet, therefore 
the leaving of vortices to infinity is possible, that is excluded for condition 



5 Three vortices on sphere 

The problem of motion of vortices on sphere is more interesting, but less investigated. 
This problem has more important applied meaning in comparison with a problem 
of motion of vortices on a plane, because it is directly connected with physics of 
atmosphere. Many substantial problems (collapse and evasion of vortices, existence 
of steady stationary (static) configurations) can be investigated on model of an ideal 
liquid. Their solutions can become a basis for more complex models (connected, 
for example, with presence of small viscosity), that describe real dynamic vortical 
processes (motion of cyclones and etc.) in an atmosphere of the Earth. 

Dynamics of two vortices on sphere is quite similar to a flat case. Here the 
general situation is the rotation around of an axis, that is passing through centre 
of spheres and analogue of vorticity centre (point located on chord, connecting two 
vorteces with radius-vector 

^_ + r 2 r 2 
r ~ r 1 + r 2 ' 

where fi,f 2 — radius- vector of the first and second vortices). As well as in a flat 
case distance between vortices remains constant. Angular speed of such rotation 
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where R — radius of sphere, M — square of distance between the vortices (square 
of the length of a chord, that connect the location of two vortices on sphere). Under 
condition of Ti = —T 2 , two vortices are moving on sphere on two identical parallels, 
located on different parties from equator, that was noticed Gromeka [3] . 

The problem of motion of three vortices on sphere was considered in work in 
that work probably for the first time the integrability of the problem was specified. 
The equations of motion of three vortices on sphere in variables of mutual distances 
(chords) and volumes can be written down ^H] as the Hamiltonian equations with 
nonlinear algebra of Poisson brackets 



4a*. A 

{Mi,Mj} — ^— , 

{M h A} = R( aj - a k )M t + R( aj + a fc )(M i - M k ) + ^(a k M k - ajMj), 

ZK 



(30) 



that also have two central functions, one of them (linear) is the same as Q, and 
the second (cubic) function has the form 



(2A) 2 + lf{M.M: + M 2 M 2 + M 3 M 3 ) - 2R 2 {M l M 2 + 
M 1 M 3 + M 2 M 3 ) + M 1 M 2 M 3 = 0, 



(31) 



and Hamiltonian (J2J). 

Let's notice, that linear approximation of structure (J5U|) coincides with the Pois- 
son bracket of three vortices on a plane. Besides, it is simple to show, that after 
replacement 

dr A 

~dt ~ 2irM 1 M 2 M 3 

the equations of a motion of vortices turn in Lotky-Volterra system of the form 

—M^ViMiiMj-Mk), (32) 



similarly to the flat case [18]. This trajectory isomorphism is also piecewise 

Geometrical interpretation for a plane, presented on Fig. 1, can be also trans- 
ferred to sphere. Curious and earlier, probably, not marked fact is that phase trajec- 
tories in variable Mi, M 2 , M 3 for a case of sphere and plane, with given intensities, 
coincide with phase trajectories of the same Lotky-Volterra system ([30)1 . Thus main 
effects in dynamics of vortices are determined by the part of phase trajectories of 
Lotky-Volterra system gets that in area A 2 (M 1? M 2 , M 3 ) > 0. The motion of vortices 
occurs only in this part, because with the approach of trajectories to border of area 
(and achievement by vortices collinear state) in the equations (|32j) it is necessary to 
change a sign of time. The form of this area is a little various for a plane and sphere 
(formula (JH) and (j3T|) ) and represented on Fig. 7. 
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Remark 3. The analogy with Lotky-Volterra system is very useful for study 
of vortices collaps problem and is related to Bolin regularization in the 11 if (the 
difference is, that in the Kepler problem is fixed a constant of energy level). Problems 
of regularization and collapse, that arise only in noncompact cases, will be discussed 
in the next part of work. 

Remark 4. The issue of reduction of three vortices problem on sphere to canoni- 
cal Hamiltonian system with one degree of freedom is complex and, probably, it could 
not be solved constructively. If one use the initial form of record llttf . then one ob- 
tain the six- dimensional Hamiltonian system, that have as the integrals Hamiltonian 
function and noncommutative set of integrals of moment, each of them is nonlinear 
function of phase variable (in difference from a case of a plane considered in \21\j). 

The canonization of the reduced system also problematic, in variable (M, A) the 
canonization is equivalent to introduction of Darboux coordinates on two-dimensional 
symplectic sheet determined by a common level of Casimir functions. It require the 
the solution of nonlinear system of the partial differential equations. 

Really, the nonlinear algebra of Poisson brackets A3U\) cannot be investigated in 
same detail, as in the flat case. Linear approximation of this structure is capable to 
make some qualitative conclusions only for a situation close to simultaneous collaps 
(i.e. when distance between vortices is small in relation to radius of curvature). 
From it, nevertheless, it is possible to make a conclusion that necessary conditions 
of simultaneous collaps of vortices on a plane are also fair for the case of sphere, 
since the influence of nonlinear terms is small enough near the collaps. 

The most simple form the bracket < f3^) has in case of three equal intensities at a 
zero level of linear Casimir function (D = 0). After transition (for linearization of 
a bracket ^EW) to canonical basis |3J) it is possible to present bracket 

{ei,e 2 } = e 3 + ~{e 2 3 -el), 

{ea.es} = ei, (33) 
{e 3 ,ei} = e 2 - e 2 e 3 . 

It is curious to note, that nonlinear terms, included in a bracket, |J are similar 
to the terms, originating in Kovalevsky integral of the Euler-Poisson equations 

One of main sections of sheet for a case T\ = r 2 = r 3 are shown on Fig. 5 in 
dependence from curvature of sphere k = 1/R. With k = 0, the symplectic sheet, that 
is determined by algebra So(3), is also sphere. The change k results in a surface, 
whose coherent component is homeomorphic to sphere. It would be interesting to 
investigate the topology of a symplectic sheet with various values of intensities and 
values of integral of the moment D, that is the limited function on sphere (on which 
vortices move). Last remark allows to reject uncompact components of symplectic 
sheets, since a motion on sphere is always finite. 

Remark 5. The general principle of the classification of vortices motions on 



11 



sphere given below is the continuation on parameter of complete moment of station- 
ary configurations that is known near D = 0. In last case the influence of curvature 
is not significant, that corresponds already investigated flat problem (it is also equiv- 
alent to consideration of linear approximation of structure \3U\) ). 

The conditions for Thomson's of configurations on sphere are similar to condi- 
tions on plane 

Mi = M 2 = M 3 . 

Conditions for existence of collinear configurations (A = A = 0) in case of sphere 
are reduced to system of three algebraic equations. Introducing constants 

«i = r 3 +r 2 , a 2 = Ti+r 3 , a 3 = r 2 +r 1; f3 x = r 3 -r 2 , p 2 = ri-r 3 , p 3 = r 2 -ri, 

we write down system of the equations as 

2(M 1 M 2 + M X M 3 + M 2 M 3 ) - (M 2 + M\ + M 3 2 ) - M ^2 M z = Q ^ 
a x M x + a 2 M 2 + a 3 M 3 - D = 0, 

(M 3 - M 2 ) (M x - M 3 ) (M 2 - M x ) (34) 
2(Ql M X + ° 2 M 2 + ° 3 M 3 + 

(/3 1 M 1 +/3 2 M 2 + /5 3 M 3 )-^ = 0, 

where D is the value of integral of the moment. The given system can be solved 
numerically. Results of numerical research of bifurcational diagrams for cases of 
various values of intensities are presented on Fig. 6. 

Remark 6. As well as in case of a plane we are interesting only in positive roots 
of the equations \3$ . Besides we should notice, that range of values of functions, 
determined by integrals of energy and moment are limited by values D m and E m 

n . D2 ( a i a 2 + a 2 a 3 + aia 3 ) 2 



CL x CL 2 Ct 3 



E m = (4R^ J] ( Z a \ + \ +a3 . 2 ) 
V- \a x + a 2 + a 3 + a x a 2 a 3 /a z k J 

In the latter case appropriate value of the complete moment is calculated by the 
formula 



d E = (4:R 2 )(a x + a 2 + a 3 ) ^ 



. y a x + a 2 + a 3 + a x a 2 a 3 /al) 



As initial approach for the solutions of the equations (jSljl with small values of 
the complete moment we chose the roots of the equations ([23)1 . Continuation of 
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bifurcational curves with increase of the moment (mutual distances) was made by 
predictor-corrector method with use of the Newton iterative procedure. 

For each set of intensities we give the bifurcational diagrams and plots of relative 
and absolute motions, on which we show the value of angular speed, corner of an 
inclination of a plane of vortices under the relation to an axis of rotation, square of 
eigenvalue, that is determining the stability of linearized system. 

Interesting effect of a spherical motion of vorteces, absent in flat case, is the birth 
of new (and in a considered case steady) collinear configurations from a problem of 
two vortices by values di = 4R 2 (cij + a*,). There is a disintegration of one vortex of 
total intensity Ti+Tj to two vortices with intensities Tj and r,- with increase of D. As 
shown in Fig. 6, these configurations aspire to merge with collinear configurations, 
turning out with continuation on parameter from a flat problem (with small D/R 2 ), 
and then cease to exist, with further increase of D. 

Let's consider bifurcational diagram for a case ai 7^ a 2 7^ 0,3 7^ a\ (see Fig. 6a). 
The Thomson's configuration (top curve) exist only by energy E < Et- With 
D = dx, by reaching maximum possible value E T < E m , Tomson's configura- 
tion merges (point A) with collinear configuration, determined by the most bottom 
branch (appropriate to bottom branch on Fig. 2 bifurcational diagram for plane). 
Thus all three vortices lay in an equatorial plane, but distance between them are 
not equal each other. Passing through a maximum of energy E m by D = ds, this 
configuration evolve futher with increases of D to a problem of two vortices (with 
D = d 3 ). Two other collinear configurations appropriate to collinear curves on Fig. 2 
for the plane, merge with collinear solutions, arising from problems of two vortices 
with values of the moment di,d 2 (disintegration of one vortex). Such arising con- 
figurations are absent in case of a plane. The occurrence could be understand from 
Fig. 8, that shown take off physical area from borders, on that the moment achieved 
value M k = 0. 

In if the ratio a\ = a 2 7^ a 3 is true, there is a merge not only Thomson's 

and collinear branch in a point A, but also collinear configurations, arising by D = 
d\ = d 2 (see. Fig. 6b). 

At last for three equal intensities all possible configurations merge together in a 
point A by the maximal values energy E = E rn = Et and moment D = dr- 

In an absolute motion with increase of the moment (mutual distances) angular 
speed of rotation of Thomson's configurations monotonously decreases (see Fig. 10), 
and corner of an inclination normal of a triangle to an axis of rotation grows 
monotonously also from zero (as for a plane) up to value 7r/2 by the moment of 
merge of Thomson's and collinear configurations (Fig. 9), by exception of a case 
equal intensities, when an inclination of an axis change not. 

Remark 7. The only possible motions of stationary configurations (i.e. con- 
figurations, with distance between vortices do not varies) are the rotations around 
of some motionless axis. It is a consequence of representability of the equations of 
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vortical dynamics as the first order equations relatively to dynamic variables. 

In contrast of Thomson's configuration, axis of rotation of collinear configura- 
tions always lays in a plane of vortices and does not vary by change values of the 
moment. 

With increasing of the moment the diagrams of angular speed are slumped 
monotonously and disappearance with merging. It is necessary to note large size 
of angular speeds of vortices, arising at the moment of birth of a new vortex from 
a problem two vortices. Within the framework of the accepted model this increase 
of angular speed concerns to different trajectories, but if is present "weak" dissipa- 
tion, that it's possible observe slow evolution both values of energy and moment. 
Certainly, there are the additional condition existing of such effects in a physical 
situation more important is the stability of appropriate stationary motions. 

Part of collinear configurations, occurring from similar configurations on a plane, 
are unstable already in linear approach, as it is visible from Fig. 11, because they 
being analogues configurations on a plane. However collinear configurations, ap- 
pearing from a problem two vortices, are steady. The nature of this stability is 
well visible from geometrical interpretation presented on Fig. 8. It is possible, that 
the phenomena of such grades occurring in an atmosphere of the Earth (that have 
obviously dissipation). That is crucial for occurrence of various catastrophic pro- 
cesses (type hurricanes), accompanying sharp reorganization of dynamics of vortical 
formations. Inverse process of collaps (merge) of two vortices, that is impossible for 
model of an ideal liquid, in case of small dissipation and reduction D can result in 
formation of atmospheric vortices with large angular speed of rotation. 

The Thomson's solutions are steady up to the moment of passage through static 
configuration. The analogue of the formula (J25J) for Thomson's configurations on 
sphere is the expression 



That shows, that Tomson's configuration is unstable with value the moment 



appropriate to the maximal value, for such configurations. 

For of one negative intensity (Fx < 0, — Ti > T 2 + T^) bifurcational 

diagram is given on Fig. 6d. In this case behaviour bifurcational curves by increase 
D is similar to behaviour by already considered situations. Existing in case of a plane 
Thomson's and collinear configurations merge in a point A (see Fig. 6d), and after 
this disappear by merging from one of collinear branches, that birth from a problem 
t of two vortices. Difference from a case only positive by increase D intensities 
is displayed in existence of collinear solutions, not limited on energy from above. 
These solutions occur also due to stationary configurations of a problem of two 




(35) 



D > 3R 2 ( ai + a 2 + a 3 ), 
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vortices. The change of parameters of an absolute motion of vortices is qualitative 
by nothing differs from a case positive intensity except for that, that all of collinear 
configuration are steady, in that time as Thomson's — is unstable (Fig. 11). 

In summary we shall explicitly allocate the basic differences of a spherical case 
from flat, arising with increase of the complete moment D : 

1. arising of Thomson's and collinear configurations on sphere; 

2. arising of steady collinear configurations from a problem of two vortices rotat- 
ing with large (infinite) angular speed in the moment of occurrence; 

3. an inclination and evolution of a plane Thomson's configurations; 

4. existence of static configurations on sphere. 

Researches of other opportunities of a motion of three vortices on a plane and to 
sphere resulting to collaps and which is running up to trajectories will be are given 
in the following part of work. 

The authors thank I.S.Mamaev and N.N.Simakov for useful discussions and help 
in work. The work is carried out under the support of Russian Fond of Fundamental 
Research (96-01-00747) and Federal programm "States Support of Integration High 
Education and Fundamental Science" (project No. 294). 
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Figure 1: Dark part corresponds to nonphysical area on a plane, that is set by linear 
integral D = const and limited by conditions Mk > 0, for a compact case: a) All of 
intensities are positive; b) One of intensities is negative. The point A corresponds 
to the Tomson's solution; F — collilear solutions; B,C,D — points, in with which 
two vortices are merged in one. 




Figure 2: The view of bifurcation diagram for three vortices on a plane in compact 
case: a) With positive intensities; b) One negative intensity. The dot-dashed line 
corresponds to the Thomson's solution, continuous lines — to collinear configura- 
tions. 
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Figure 3: Dependence of a squares of eigenvalues, determining stability in linear ap- 
proach from integral of complete moment for stationary configurations in a compact 
case: a) With positive intensities; b) One negative intensity. The dot-dashed line 
corresponds to the Thomson's solutions, continuous lines — to collinear configura- 
tions. 




Figure 4: Phase plane for cases: a) Positive equal intensities; b) Case a x = a 2 ^ a 3 
c) Three positive various intensities; d) One negative intensity. 
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Figure 5: Deformation of section of symplectic sheet for linearized algebra of vortices 
on sphere depending on change of sphere's radius of curvature. 




Figure 6: The bifurcational diagram for cases: a) Various positive intensities; b) 
Positive intensities (two coincide); c) Equal positive intensities; d) One negative 
intensity. The designations of curves coincide with Fig. 2. 
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Figure 7: Geometrical picture is isomorfic with a Lotky-Volterra problem. A part 
of a phase plane of Lotky-Volterra, that is limited by a contour, corresponds to 
physical area of a problem of three vortices: a) on a plane; b) on sphere with a set 
of identical intensities. 
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Figure 8: Geometrical picture of birth of collinear solutions from a problem of 
two vortices. Consecutive figures a), b), c) show change of points of a contact of 
borders of physically of allowable area with increase of the complete moment. The 
designations correspond Fig. 1. 




Figure 9: Angular speeds for cases: a) Various positive intensities; b) Positive in- 
tensities (two coincide); c) Equal positive intensities; b) One negative intensity. 
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Figure 10: Corners of an inclination for cases: a) Various positive intensities; b) 
One negative intensity. 




Figure 11: Factor of stability A 2 for cases: a) Various positive intensities; b) Positive 
intensities (two coincide); c) Equal positive intensities; d) One negative intensity. 
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